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ABSTRACT 

Recently, the CANGAROO-II telescope detected diffuse TeV gamma-ray 
emission from a nearby edge-on starburst galaxy, NGC 253. The emission mecha- 
nism is discussed in this report. We review the emissions of radio-to- TeV gamma- 
rays from NGC 253, and present a model of the non-thermal emissions due to 
synchrotron radiations and inverse Compton scatterings. A halo model success- 
fully explains the multiband spectrum of NGC 253. 

Subject headings: cosmic rays — galaxies: halos — galaxies: individual(NGC 253)— 
galaxies: starburst — gamma rays: theory 



1. Introduction 



The CANGAROO-II imaging atmospheric Cherenkov telescope has detected TeV gamma- 
ray emission from a nearby edge-on starburst galaxy, NGC 253 (Itoh et al. 2002a,b). This 
TeV emission is spatially extended and temporally steady, whose nature is different from 
that of previously observed extragalactic objects of the active galactic nuclei (AGN) class. 
This is the first detection of TeV gamma-rays from a normal-sized spiral galaxy like our 
Galaxy. We can learn how high-energy particles propagate on the galactic scale from the 
observations of NGC 253 by TeV gamma-rays. 

The acceleration and propagation of galactic cosmic rays (GCRs), which are observed 
directly near the earth, are among the big topics in physics and astrophysics. The distribution 
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of GCR with an energy of less than ~10 GeV in the Galaxy has been surveyed by tracing 
diffuse gamma-rays which are produced by the interaction of GCR with the interstellar 
medium(Hunter et al. f997). The resultant distribution of diffuse gamma-rays as a function 
of the galactic longitudes and latitudes has been argued based on the scenario of the SNR 
origin of GCR. In the TeV energy range (LeBohec, S. et al. 2000; Aharonian et al. 2001), 
searches for diffuse gamma-rays in our Galaxy have been also carried out to obtain the 
spatial distribution of TeV cosmic rays. It is, however, difficult due to the limitations of 
the present techniques. In some respects, studies of nearby galaxies which are near enough 
to obtain detectable fluxes from outside are more reasonable than that of our Galaxy from 
inside. A further question naturally arises as to whether other galaxies also harbour cosmic- 
rays like ours. Photons from radio to gamma-rays produced by cosmic-rays would reveal 
their existence. The electrons of GeV energies associated with edge-on galaxies have been 
known by radio observations (Berzinskii et al. 1990; Beck 1997). EGRET demonstrated clear 
evidence for the existence of GeV particles in LMC (Sreekumar et al. f992); however, no 
evidence for such high-energy particles as TeV in normal galaxies had been known, except 
for that from our Galaxy. Obviously, investigating nearby galaxies like NGC 253 is a key. 

In Section 2 we review the emissions of radio-to- TeV gamma-rays of this galaxy. In 
Section 3 we present a model which explains the observed multiband spectrum and the 
spatial distribution of cosmic rays in NGC 253. We discuss the propagation of high-energy 
particles into the halo, the reacceleration of particles due to the galactic wind, and the total 
cosmic-ray energy of NGC 253 in Section 4. 

2. Multi- wavelength Spectrum 

The distance to NGC 253 is as close as 2.5 Mpc (Vaucouleurs 1978) with a visual size 
of 0.3° from Earth. This is one of the best targets where the cosmic-rays density is expected 
to be high due to its high supernova rate (Volk et al. f996). Here, we summarize previous 
studies at NGC 253 in various photon energies. The compiled fluxes and upper limits on 
NGC 253 are shown in Fig. I. 

The EGRET instrument on board the Compton Gamma Ray Observatory(CGRO) has 
reported stringent upper limits for GeV gamma-rays (Sreekumar et al. f994; Blom et al. 
1999). On the other hand, the OSSE on board CGRO has claimed the detection of 50^200 
keV gamma-rays (Bhattacharya et al. 1994). The PDS on board BeppoSAX gives a contra- 
dicting upper limit on the hard X-ray emission just below OSSE's energy band (Cappi et al. 
1999). 
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At TeV energies, CANGAROO-II has detected diffuse emission. The size of the emitting 
region is similar to, or larger than, the optical image of NGC 253 (Itoh et al. 2002a,b). A 
simple power-law fit gave a soft spectral photon index of 3.74 ± 0.27(Itoh et al. 2002b), 
suggesting a cutoff. 

The radio continuum emission consists of three main components: a central region, a 
disk component, and a halo extending ~ 9 kpc (Hummel et al. 1984; Carilli et al. 1992). 
The radio spectral energy index is —0.7 ± 0.1 in the disk, and it is steepened to ~ —1 in 
the halo. The fractional polarization increases with the distance from the plane (less than 
4% in the disk, up to 15% in the halo in 1.4GHz (Carilli et al. 1992), and reaching 40% 
at 10GHz (Beck et al. 1994)). The turbulent magnetic fields estimated from the fractional 
polarization are ~ 17fiG in the disk and ~ QfiG in the halo (Beck et al. 1994). The radio 
emission in the central region of ~ 150 x 70pc shows a flat spectrum with an energy index 
of ~ —0.43 (Hummel et al. 1984). From high-resolution radio observations of the inner 200 
pc of NGC 253, Ulvestad & Antonucci (1997) detected ~ 64 compact radio sources, half of 
which are in the HII regions with flat spectra. The other half are presumably synchrotron 
supernova remnants. 

The X-ray structure of NGC 253 is complex, and consists of point sources, hot inter- 
steller matter (ISM) in the disk and the central region, and the halo. The diffuse X-ray 
emission from the halo was first detected with Einstein, and followed up by ROSAT, ASCA, 
BeppoSAX, XMM-Newton, and Chandra (Watson et al. 1984; Fabbiano 1988; Pietsch et al. 
2000; Ptak et al. 1997; Pietsch et al. 2001; Strickland et al. 2000, 2002; Cappi et al. 1999; 
Weaver et al. 2002). X-ray spectral studies of the disk and the halo show the existence of 
optically thin thermal plasmas at temperatures of kT = 0.1 ~ 0.4 keV. Chandra and XMM- 
Newton observations of the central region of NGC 253 reveal the existence of a heavily 
absorbed high-temperature thermal plasma component (kT ~ 6 keV) with an ionized iron K 
emission line (Pietsch et al. 2001; Weaver et al. 2002). This absorbed component is also seen 
in the spatially integrated spectra over the whole galaxy obtained with the previous X-ray 
satellites (Ptak et al. 1997; Cappi et al. 1999; Ohashi et al. 1990). Based on a magnetic 
field strength of 270 /zG in the central region, estimated from the standard minimum energy 
assumption, Weaver et al. (2002) suggested that the inverse Compton process makes only a 
minor contribution to the hard X-ray emission. 

3. Inverse Compton Emission Halo Model 

We first discuss the emission mechanism at gamma-ray energies from the multi- wavelength 
spectrum, as shown in Fig. 1. Goldshmidt & Rephaeli (1995) reported that the gamma-ray 
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emission detected by OSSE (the blank circles in Fig. 1) was attributed to the Inverse Comp- 
ton Scattering (IC) of ubiquitous far-infrared(FIR) photons from dusts observed around the 
center of the galaxy. However, a simple extrapolation to the GeV range exceeds the 2 a 
upper limit of 3.4 x 10~ 8 cm~ 2 s _1 (Blom et al. 1999) obtained by EGRET (the red trian- 
gles). The multi- wavelength spectrum suggests that the origin of TeV gamma-rays should 
be different from that of sub-MeV gamma-rays. 

If we assume that the energy spectrum of particles is a power law with an exponential 
cut, and that the power-law index is softer than —2, the index of the differential photon 
spectra in GeV energies due either to the emissions of 7r° decay produced by the collisions 
of high-energy protons with the matter or non-thermal bremsstrahlung of the high-energy 
electrons also becomes softer than —2, leading to a contradiction with the EGRET upper 
limit. 

The remaining possibility of the emission mechanism of GeV and TeV gamma-rays is 
IC scattering. When the power-law index of the electron energy spectrum is assumed to be 
~ —2 , the index of the differential photon spectrum of the IC becomes nearly —1.5. 

The very large radio halo of ~ 9 kpc (Carilli et al. 1992) would suggest the existence 
of a population of very high-energy electrons which emit our TeV gamma-rays other than 
those concentrated near to the central or disk region of the galaxy, which might emit photons 
observed by OSSE. On the other hand, there are sure target photons for IC, such as CMB as 
well as very abundant FIR photons. Here, we quantitatively examine whether these electrons 
and photons are responsible for TeV emissions. 

We assume that the total number spectrum of electrons at the source has the following 
form: 



where N e0 is an electron density factor, 7 the Lorentz factor, and 7m the maximum Lorentz 
factor of the electrons (a cutoff), respectively. The spectrum of the IC emission is given by 



where D is the distance to the source, n p h the number spectrum of the soft photon, c the 
speed of light, e the energy of the target photons, E 1 the energy of the scattered photons, E 
the energy of the electrons, and <r(E~ n e, E) the Klein-Nishina cross section for IC scattering, 
respectively. 

We take the FIR photon spectrum model, in which the dust absorption efficiency is 
proportional to the photon frequency, as described in Goldshmidt & Rephaeli (1995). The 
total FIR luminosity is L IR = 1.97 x 1O 1O L , with 1.45 x 1O 1O L from cool dust having a 



^(7) = V e0 7" p exp(-7/ 7m ), 



(1) 
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temperature T c = 36.8 K and 0.52 x 1O 1O L from warm dust with T w = 172.8 K, respectively. 
Here, we assume that the photon field is isotropic. The energy density of the photon field 
C/tr at radius R s is approximated as £/ir = Li R /(nR 2 c). We simply adopt the average value 
of £7ir in a spherical region with a size of R s as the photon density. The average photon 
density (< U m >) is f* s U m (r)47ir 2 dr /(4ttR s 3 /3) = 3U m (R s ). The optical photon field 
is assumed to be a diluted blackbody with a temperature of 6000 K and a luminosity of 
L B = 2.65 x 1O 1O L (Pence 1980). Assuming the value of the electron spectrum index p, we 
could determine the two parameters iV e0 /(47i\D 2 ) and 7 m for the electron spectrum as fitting 
the CANGAROO-II flux and keeping the EGRET upper limit. 

Here, we adopt a one- zone model with two different sizes for simplicity. First, we adopt 
R s = 10 kpc, suggested by the size of TeV emitting region (Itoh et al. 2002a), with assuming a 
cosmic-ray halo. We considered three populations of photons as the target of the IC process, 
IR, optical, and CMB. The energy densities for each population are < [/"ir >= 1.6 eV cm" 3 , 
< U opt >= 2.2 eV cm" 3 , and U C mb = 0.26 eV cm" 3 , respectively. Although < U opt > is 
greater than < [/tr >, its contribution was suppressed by the Klein-Nishina effect. In Fig. 
1, the solid and dashed lines represent the synchrotron and IC emission models with the 
power-law index of the electrons being p = 2 and 2.2, respectively. Here, the value of p was 
chosen based on constraints by the shock acceleration (p > 2) and reconciliation between the 
EGRET upper limits and the CANGAROO-II TeV results (p = 2.4) (Itoh et al. 2002a,b). 
The fluxes of the TeV gamma-rays and the upper limit of EGRET restricted the parameter 
choices of the electron spectrum: N e0 / '(An D 2 ) = 6.8 x 10 8 cm" 2 and 7 m = 2.5 x 10 6 (1.3TeV) 
for p = 2, N e0 /(4nD 2 ) = 8.4 x 10 9 cm" 2 and 7 m = 3.6 x 10 6 (1.9TeV) for p = 2.2, and 
N e0 /(4nD 2 ) = 8.5 x 10 10 cm" 2 and 7 m = 4.7 x 10 6 (2.4TeV) for p = 2.4, respectively. The IC 
fluxes at ~ 2 keV for p = 2.2 ~ 2.4 are predicted to be the same as the Chandra X-ray data 
in the halo region. This X-ray data imply these halo models favor the value of p < 2.2 ~ 2.4. 
We find the total energies and average energy densities of the electrons to be within R s = 10 
kpc: 5.9 x 10 54 erg and 0.03 eV cm " 3 for p = 2, 2.4 x 10 55 erg and 0.12 eV cm " 3 for 
p = 2.2, and 1.3 x 10 56 erg and 0.66 eV cm ~ 3 for p = 2.4, respectively. The strength of 
the magnetic field (B) were estimated to be 2.5/xG for p = 2, 1.7/xG for p = 2.2, and 1.3;uG 
for p = 2.4, respectively, by adjusting the calculated synchrotron emissions to the radio halo 
data. These values are smaller than the estimations obtained by the fractional polarization 
in radio measurements (Beck et al. 1994). 

On the other hand, we examined the case for a disk component of cosmic-rays. We 
adopted a higher average photon density within R s = 3 kpc. The values of the average 
photon energy densities are < Uir >= 18 eV cm" 3 and < U opt >= 24 eV cm" 3 . In Fig. 1, 
the dotted lines represent the emission models from a disk with the power-law index of the 
electrons of p = 2.4, which was estimated based on the radio data of the disk. The parameters 
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of the electron spectrum are N e0 /(4irD 2 ) = 7.7 x 10 9 cm- 2 and 7 m = 4.7 x 10 6 (2.4 TeV). We 
found the total energy and average energy density of the electrons within R s = 3 kpc to be: 
1.2 x 10 55 erg and 2.2 eV cm ~ 3 . The strength of the magnetic field, B = 8.2/xG, reproduces 
the radio fluxes from the disk. In the case of a soft electron spectrum with p = 2.4, fittings 
for the CANGAROO-II fluxes and the upper limits of EGRET are rather difficult, as shown 
in Fig. 1. The high-energy electron halo model explained better the radio and GeV-TeV 
emission. 

The halo model, however, cannot describe any emissions in the keV-MeV energy region. 
Goldshmidt & Rephaeli (1995) introduced the IC model with a very localized FIR and high- 
energy electrons in order to explain the OSSE result. Their IC model should require a cutoff 
in the energy spectrum of electrons to reconcile with the EGRET upper limits. We could 
reproduce the OSSE flux by adopting an average photon density of < Um >= 1.8 x 10 3 eV 
cm' 3 within R s = 0.3 kpc, and taking the following parameters of the electron spectrum: 
for example, in case of p = 2.3, N e0 /(AivD 2 ) = 4.5 x 10 9 cm" 2 and 7m = 1.2 x 10 4 (6.2 GeV). 
The total energy and average energy density of the electrons are 8.5 x 10 54 erg and 1.6 x 10 3 
eV cm ~ 3 , respectively. Such a high density and low cutoff energy may be reasonable, 
considering the starburst characteristics and FIR density around the center of this galaxy. 
This is also consistent with the EGRET upper limits, as shown in Fig. 1 (the long-dashed 
line). The synchrotron emission was also reproduced along with the strengths of the magnetic 
field, B = 9.6//G, which is almost the same value as that adopted by Goldshmidt & Rephaeli 
(1995). However, the model IC flux at or above ~ 2 keV is more than one order of magnitude 
higher than the Chandra and XMM X-ray data, which were corrected for the photoelectric 
absorption in both the Galaxy and NGC 253. The model IC flux in hard X-ray energy band 
is also higher than the BeppoSAX PDS X-ray data, while these data were not corrected 
for the absorption which was estimated to be negligibly small in the higher energy. We 
concluded that some modifications are necessary for this model. 

Hitherto, we considered only the possibility of a leptonic origin of the observed TeV 
gamma rays; however, we should check before reaching any final conclusion whether hadronic 
models can explain CANGAROO-II's positive detection consistently with the EGRET upper 
limits. In this estimation of gamma-ray emission by 7r°-decay(Mori 1997), we assume the 
proton-to-electron ratio to be 300, which is rather high compared to that of GCRs (~ 100) 
and accordingly is a favourable case for hadronic models. We assume the parameters for 
electrons to be p = 2 with a higher cutoff energy of 10 TeV and N e0 /4irD 2 = 6.8 x 10 8 cm~ 2 , 
which we have shown previously. We also assumed the total mass of the gas which interacts 
with protons to be 3 x 1O 9 M (Puche et al. 1991; Houghton et al. 1997). The expected 
gamma-ray flux is shown in Fig. 1 as the dot-dashed line. It can be clearly seen that 
hadronic models fail to explain the observations, even in the favourable case. 
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4. Propagation and Acceleration 

The high-energy electrons in the halo of the galaxy explain most of multiband spectrum 
of NGC 253 quite well. Here, we discuss the origin of these cosmic-ray electrons and the 
total cosmic-ray energy, which should be compared with our Galaxy. 

The biggest question is whether high-energy electrons can propagate out to a distance of 
~ lOkpc from the disk of the galaxy without any severe energy losses. Electrons may diffuse 
out to a distance of Rl ~ 2(/tt cooi ) 1 / 2 within their cooling time (t coo i) due to synchrotron and 
IC losses, where t coo i is inversely proportional to a sum of the energy densities of the magnetic 
field and the photons in ambient space, and is also inversely proportional to the energy of the 
electrons. When the diffusion coefficient of electrons (k) of order 3 x 10 29 (i?/GeV) - 6 cm 2 s _1 
is used, we found that electrons of 1 TeV which emit sub-TeV gamma-rays can propagate to 
a distance of ~ 9 kpc, assuming the strength of magnetic field to be 2.5/xG, and the energy 
density of FIR photon to be C/ir(10 kpc). The assumption regarding the energy dependence 
of k is crucial for this estimation. 

The diffuse X-ray emission extending from the nuclear region toward the halo is thought 
to be due to the galactic wind (super-wind), powered by supernovae occurring in the nuclear 
starburst region (Chevalier & Clegg 1985), which has been successfully described by many 
papers based on hydrodynamical simulations (Tomisaka et al. 1988; Strickland & Stevens 
2000; Suchkov et al. 1994). Continuous supernova explosions heat up the ISM in the nuclear 
starburst region, and creates a bubble of high-temperature plasmas with T ~ 10 8 K. The 
bubble expands while sweeping and pushing cool ISM and ambient gas up towards the halo. 
A galactic wind running freely along the galactic minor axis is formed around 10 7 yr after 
the start of the starburst activity. Wind with a velocity of 2000 ~ 3000km sec -1 finally hits 
the ambient cool gas and forms a shock in the halo. This picture explains that the image of 
the synchrotron radio halo at 0.33GHz generally matches the diffuse X-ray emission (Pietsch 
et al. 2000). This shock in the halo may accelerate particles by the Fermi process. 

The maximum energy of electrons accelerated by this process is estimated to be 8 ~ 
25 TeV if we assume a shock velocity of ~ 2000km sec -1 , a magnetic field of 2 fiG, and a 
starburst age of ~ 10 7 yr by equating the acceleration rate with the cooling rate due to 
the synchrotron and IC losses. The estimated maximum energy of electrons agrees well with 
that determined from TeV gamma-ray spectrum by CANGAROO observation. Another 
possibility of particle acceleration in the termination shock of the galactic wind has been 
proposed by Jokipii & Morfill (1985). These two cosmic-ray acceleration processes may be 
related to the rather hard spectrum required for the electrons in the halo. 

Finally, we discuss the total cosmic-ray energy in the halo. The electron energy density 
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was obtained to be 0.03 ~ 0.12 eV cm -3 , which is one order higher than that of our Galaxy. 
If we assume the proton to electron ratio to be 100, this implies that the total energy of 
cosmic rays in NGC 253 may amount to 5.9 x 10 56 ~ 2.4 x 10 57 ergs, which is one hundred- 
times higher than that for our Galaxy. Note that the value for our Galaxy was derived 
assuming a disk population of cosmic-rays. If our Galaxy has a halo with the same cosmic- 
ray density, the deviation becomes smaller. The starburst age and the supernova rate in 
NGC 253 were estimated to be ~ 10 7 yr and ~ OSyr^ 1 , respectively (Rieke et al. 1988). 
Thus, taking the ambiguities (e.g. the proton to electron ratio) into account, the starburst 
activity characterizing NGC 253 may explain the large amount of total energy of cosmic 
rays. 

We gratefully acknowledge the CANGAROO-II team for discussions and suggestions, 
and thank Dr. M. Mori for his support in providing us with his code for calculating the 
gamma-ray spectrum from the pion decay. We also thank Prof. W.Pietsch for his kindly 
providing the information on the XMM spectrum. 
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Fig. 1. — Multi-band spectrum of NGC 253. The black squares were obtained by 
CANGAROO-II. The X-ray data were corrected for photo-absorption in both the Galaxy 
and NGC 253, except that no correction was applied to the BeppoSAX hard X-ray data. 
We note that the BeppoSAX, OSSE, and EGRET data were not able to spatially resolve 
NGC 253, in contrast to the other data. The lines shown were obtained by estimations 
described in the text. 



